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Abstract 

Alkali halide (100) crystal surfaces are poorly wetted by their own melt at the triple 
point. We carried out simulations for NaCl(lOO) within the well tested BMHFT 
model potential. Calculations of the solid-vapor, solid-liquid and liquid-vapor free 
energies showed that solid NaCl(lOO) is a non-melting surface, and explain its bad 
wetting in detail. The extreme stability of NaCl(lOO) is ideal for a study of the 
nanofriction in the high temperature regime, close to and even above the bulk 
melting temperature (Tm). Our simulations reveal in this regime two distinct and 
opposite phenomena for plowing and for grazing friction. We found a frictional 
drop close to Tm for deep ploughing and wear, but on the contrary a frictional rise 
for grazing, wearless sliding. For both phenomena we obtain a fresh microscopic 
understanding, relating the former to "skating" through a local liquid cloud, the 
latter to softening of the free substrate surface. It is argued that both phenomena, 
to be pursued experimentally, should be much more general than the specific NaCl 
surface case. Most metals in particular possessing one or more close packed non- 
melting surface, such as Pb, Al or Au(lll), should behave quite similarly. 
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1 Introduction 



The high temperature surface properties of alkali halide crystals, particularly 
of NaCl, are very unusual as these solids are incompletely wetted by their own 
melt. [1] We investigated theoretically these surfaces and their uncommon be- 
havior, with a multiple goal. The first goal was to uncover the physical reasons 
for the poor wetting. The second was to exploit the availability of the standard 
empirical potential parametrized long ago by Tosi and Fumi (BMHFT) for a 
first quantitative characterization of all the interfaces - and particularly their 
structure and their temperature dependent interface free energies - in a well 
defined solid-liquid-vapor system. Achieving this kind of goal is of particular 
interest in the general context of high temperature capillarity, a field where 
such a detailed understanding is usually unavailable. A third goal was to ver- 
ify whether the solid NaCl (100) surface exhibits surface non-melting, as the 
partial self-wetting would imply on thermodynamic grounds. A fourth goal 
was to study the NaCl liquid surface, extracting from its temperature depen- 
dent surface tension the surface entropy, and comparing it. e.g., with the solid 
surface entropy at the melting point. Again, it seems that this comparison is 
generally unavailable. A fifth and more practical type of goal was finally to 
put to use the understanding obtained for the non-melting solid surface, in 
the field of nanofriction. In particular, our aim is to provide a first theory and 
simulation approach to both ploughing and grazing friction of hard nanotips 
on a solid surface near melting. 

In the following we summarize the results of work that appeared in Refs. [2H31H1I5] . 



2 Partial wetting of NaCl by its own melt 

Our classical molecular dynamics simulations with the Tosi-Fumi potentials 
confirm that in this model crystalline NaCl(lOO) is only partly wetted by a 
(nano) droplet of molten NaCl. [2J Moreover, as expected, solid NaCl (100) is 
confirmed to be a non-melting surface, stable without any precursor signals of 
melting up to the bulk melting point. In a metastable state, and in the absence 
of defects, it can even be overheated by as much as 100K above the melting 
temperature (Fig. [T]). Extracting from the simulations the properties of the 
three interfaces - solid-vapor, liquid-vapor, and solid-liquid - we calculated 
and analyzed their free energies (Fig. [2]) and found that the surface non- 
melting of NaCl(lOO) and the resulting partial self-wetting stem from three 
separate reasons: 

(i) Solid NaCl(lOO) is an exceptionally stable solid surface (Fig. [Q, with a 
strongly decreasing surface free energy at high temperature (Fig. Sta- 
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bility in this regime is enhanced by extremely large anharmonicities. At the 
melting point T M , we found 7sv = 103 mJ/m 2 . 

(ii) The solid-liquid interface is spatially sharp (Fig. [3]), and is energetically very 
expensive (7sl — 36 mJ/m 2 ). This can be attributed to large structural 
differences between solid and liquid, in particular a 27% density difference 
between the two. 

(iii) The liquid-vapor interface free energy (surface tension) is relatively high 
(7l V — 102 mJ/m 2 ), actually almost identical to the solid- vapor free energy 
at the melting point. The reason for this high surface tension is traced to 
a surprising deficit of liquid surface entropy, which we calculate to be a 
factor ~2.5 lower than that of the solid surface. Low entropy implies some 
kind of short range order at the molten salt surface. We highlighted and 
characterized in particular an incipient pairwise molecular charge order in 
the outermost regions of the molten salt surface (Fig. 0J. [3"P] 



3 High temperature nanofriction 

The unusual stability of the alkali halide surface against melting suggests 
adopting it as a natural toy system for the study of high temperature AFM 
tip-surface nanofriction studies. A non-melting surface will not automatically 
liquefy under a tip, (as would instead be the case for a regular, melting sur- 
face [6]) even at temperatures very close to bulk melting. For this reason, 
NaCl(lOO) and other alkali halide surfaces should provide good testing grounds 
for a first study of high temperature nanofriction, a potentially important area 
which is still largely unexplored. In the absence of any experimental data on 
high temperature nanofriction, we conducted exploratory molecular dynamics 
sliding friction simulations of hard (diamond) tips on NaCl(lOO), both in the 
heavy ploughing, wear-dominated regime, and in the light grazing, wearless 
regime. 

The simulated ploughing friction of a sharp penetrating tip is very large, as 
expected for scratching a hard solid (Fig. [5]). However, it shows as a function 
of temperature a strong drop near the melting point. This high temperature 
friction decrease can be assimilated to skating of the tip over the hot solid. As 
in ice skating, the tip is surrounded by a local liquid halo, moving along with 
the tip as it ploughs on. [5] 

At the opposite extreme, we found that simulated grazing friction of a hard 
flat tip over NaCl(lOO) behaves just the other way around (Fig. E]). Unlike 
ploughing, low temperature grazing friction is initially very small, on account 
of the incommensurability of the two respective lattice planes. Upon heating 
however the grazing friction surges just near the melting point. This rise is seen 
as the analog of the celebrated peak effect in the critical current close to the 
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upper critical field in the mixed state of type II superconductors. The dramatic 
softening of the flux lattice just before its disappearance causes a sudden 
increase of sliding friction of fluxons against the pinning impurities. [7f5] In 
our case the softening of the NaCl surface lattice, still solid at the melting 
temperature, reflects in an increase of sliding friction just before spontaneous 
breakdown of the solid surface above the melting point. [5] 

4 Conclusions 

This contrasting behavior of friction close to the melting point - skating when 
scratching, braking when grazing - should be more general than the specific 
instance of NaCl (100) where we uncovered them. Metals in particular usually 
possess at least one close packed non-melting surface [S] such as the (111) 
or the (110) surfaces of fee or bee metals respectively, which is non-melting 
and should behave similarly to NaCl(lOO). We moreover can expect possi- 
ble modifications of the high temperature nanofriction behavior by adhesive 
forces between tip and surface, as well as by application of a bias between 
(conducting) tips and surfaces. 
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Fig. 1. Snapshot of a simulation of a NaCl(lOO) slab, at 1120 K ~ T M + 50 K 
(the calculated bulk melting temperature is Tjyj = 1066 K, in good agreement 
with experimental data). The (100) surface is normal to the z direction, indicated 
in the figure. In our simulations at 1120 K, NaCl(100) remains crystalline in a 
metastable state for at least 200 ps, confirming its non-melting nature. Upon in- 
creasing temperature, crystalline NaCl(100) persists up to a surface spinodal tem- 
perature T ss ~ 1215 K ~ T M + 150 K. 





Fig. 2. (a) The solid surface free energy 73 v from thermodynamic integration (solid 
line) and from the effective harmonic approximation (dots) [JJ. (b) Liquid NaCl 
surface tension 7lv, compared to the experimental value 'y^f ■ The solid surface free 
energy 7sv is reported in the same plot. 7^ v is the surface tension re-calculated by 
artificially removing the correlations between atoms at the outer surface. |3|4j Once 
surface molecular order is removed in this way, the surface entropy rises and the 
surface tension drops. Solid NaCl(100) would be completely wet by this artificial 
liquid. 
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Fig. 3. Simulated NaCl(lOO) solid-liquid interface at the melting point. The interface 
is quite sharp. Its free energy is calculated to be 7sl = 36±5 mJ/m 2 . 




Fig. 4. Snapshot of a simulation of a liquid NaCl surface (liquid-vapor interface), at 
1250 K well above the melting point and even above the surface spinodal tempera- 
ture T ss . Note the very pronounced fluctuations in the instantaneous surface profile. 
This picture, suggestive of a low surface tension, high entropy surface, is in apparent 
contradiction with the massive non- wetting of solid NaCl(lOO). Indeed, we found in 
simulation |3)4j . an incipient molecular surface ordering which lowers the surface 
entropy by a factor ~ 2.5. The origin of this molecular surface order is charge order, 
which, already important in bulk, plays an enhanced role at the molecular liquid 
surface. Analyzing the correlations between Na + and Cl~ ions at the surface, one 
finds that the outermost layer is rich in NaCl monomers and dimers, as highlighted 
in the figure. 
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Fig. 5. Top panel: initial stage picture of ploughing friction of a hard, sharp tip 
on NaCl(100), 100 degrees below the melting point. The tip is modeled as a rigid 
diamond apex, 26 A long. The tip indentation depth is ~ 6 A. Note the fast-healing 
furrow behind the tip. Bottom panel: averaged ploughing frictional force (F x ) as a 
function of temperature. Note the friction drop at around T = Tm — 150 K, which is 
due to "skating" of the tip. 
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Fig. 6. Top panel: snapshot a grazing friction simulation. The tip is modeled as 
rigid hard diamond tip, exposing a flat facet of about 132 A 2 . The tip contacts 
the substrate by mere adhesion forces, which we estimate to be ~ 1 nN. Bottom 
panel: average frictional force of the grazing tip on NaCl(100). Note the initially 
very low friction (due to incommensurability between the two lattices), however 
sharply increasing at Tm and above. Note the contrast with the ploughing friction 
drop of Fig. [5] 
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